Superlight bipolarons and a checkerboard d-wave condensate in cuprates 
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The seminal work by Bardeen, Cooper and Schrieffer taken further by Eliashberg to the inter- 
mediate coupling solved the problem of conventional superconductors about half a century ago. 
The Frohlich and Jahn- Teller electron-phonon interactions were identified as an essential piece of 
physics in all novel superconductors. The BCS theory provides a qualitatively correct description of 
some of them like magnesium diborade and doped fullerenes (if the polaron formation is taken into 
account). However, cuprates remain a problem. Here I show that the bipolaron extension of the 
BCS theory to the strong-coupling regime could be a solution. Low-energy physics in this regime is 
that of small 'superlight' bipolarons, which are real-space mobile bosonic pairs dressed by phonons. 
The symmetry and space modulations of the order parameter are explained in the framework of the 
bipolaron theory. A d-wave Bose-Einstein condensate of bipolarons reveals itself as a checkerboard 
modulation of the hole density and of the gap below Tc. 

PACS numbers: PACS: 74.72.-h, 74.20.Mn, 74.20.Rp, 74.25.Dw 



I. INTRODUCTION: PAIRING IS INDIVIDUAL 
IN MANY CUPRATES 

Experimental!! S H S S H, SH and theoretical 
0, 0, 0, 0, Il3| evidence for an exceptionally strong 
electron-phonon interaction in high temperature super- 
conductors is now so overwhelming, that even some ad- 
vocates of nonphononic mechanisms |14| accept this fact. 
Our view is that the extension of the BCS theory towards 
the strong interaction between electrons and ion vibra- 
tions describes the phenomenon naturally. High temper- 
ature superconductivity exists in the crossover region of 
the electron-phonon interaction strength from the BCS- 
like to bipolaronic superconductivity as was predicted be- 
fore .and explored in g reater detail after the discovery 

E3, Ha E3, EE 123, Eif . 

Polarons form the Cooper pairs, if their coupling is 
not strong enough. However, real-space small bipolarons 
are also feasible due to a polaron collapse of the Fermi 
energy. There is strong evidence for the Bose-Einstein 
condensation of small bipolarons in many cuprates, in 
particular at low doping. Let us, for example, estimate a 
renormalised Fermi energy. The band structure is quasi- 
two-dimensional with a few degenerate hole pockets in 
cuprates. Applying the parabolic approximation for the 
band dispersion we obtain the renormalized Fermi energy 
as 



(1) 



where d is the interplane distance, and m* arc the den- 
sity of holes and their effective mass in each of the hole 
subbands i renormalized by the electron-phonon (and 
by any other) interaction. One can express the renor- 
malized band-structure parameters through the in-plane 
magnetic-field penetration depth at T « 0, measured ex- 



perimentally: 
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As a result, we obtain a parameter- free expression for the 
"true" Fermi energy as 



d 



=2 A 2 
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where g is the degeneracy of the spectrum, which may 
depend on doping. One expects 4 hole pockets inside the 
Brillouin zone (BZ) due to the Mott-Hubbard gap in un- 
derdoped cuprates. If the hole band minima are shifted 
with doping to BZ boundaries, all their wave vectors 
would belong to the stars with two or more prongs. The 
groups of wave vectors for these stars have only ID repre- 
sentations. It means that the spectrum will be degener- 
ate with respect to the number of prongs, i.e g ^ 2. Be- 
cause Eq.(3) does not contain any other band-structure 
parameters, the estimate of €f using this equation does 
not depend very much on the parabolic approximation 
for the band dispersion. 

Generally, the ratios n/m* in Eq.(l) and Eq.(2) are not 
necessarily the same. The 'superfiuid' density in Eq.(2) 
might be different from the total density of delocalized 
carriers in Eq.(l). However, in a translationally invari- 
ant system they must be the same • This is also true 
even in the extreme case of a pure two-dimensional su- 
perfiuid, where quantum fluctuations are important. One 
can, however, obtain a reduced value of the zero tempera- 
ture superfiuid density in the dirty limit, / <C £(0), where 
£(0) is the zero-temperature coherence length. The latter 
was measured directly in cuprates as the size of the vor- 
tex core. It is about 10 A or even less. On the contrary, 
the mean free path was found surprisingly large at low 
temperatures, / ~ 100-1000 A. Hence, it is rather prob- 
able that novel superconductors are in the clean limit, 
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TABLE I: The Fermi energy (multiplied by the degeneracy) 
of cuprates 

Compound T c (K) X H ,ab (A) d(A) ge F (meV) 
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to get an experimental T c , one has to have a strong cou- 
pling, A > fi c . However, one cannot increase A without 
accounting for the polaron collapse of the band. Even in 
the region of the applicability of the Eliashberg theory 
(i.e. at A < 0.5), the non-crossing diagrams cannot be 
treated as vertex corrections, since they are compara- 
ble to the standard terms, if luq/ef > 1. Because novel 
superconductors are in the nonadiabatic regime, interac- 
tion with phonons must be treated in the framework of 
the multi-polaron theory at any value of A. 

In many cases (Table) the renormalized Fermi energy 
is so small that pairing is certainly individual, i.e. the 
bipolaron radius is smaller than the inter-carrier distance. 
Indeed, this is the case, if 

£ F <7tA. (5) 

The bipolaron binding energy A is thought to be twice 
the so-called pseudogap experimentally measured in the 
normal state of many cuprates, A > lOOmeV, so that 
Eq.(5) is well satisfied in underdoped and even in a few 
optimally and overdoped cuprates. One should notice 
that the coherence length in the charged Bose gas has 
nothing to do with the size of the boson. It depends on 
the interparticle distance and the mean-free path, [Tof . 
and might be as large as in the BCS superconductor. 
Hence, it is incorrect to apply the ratio of the coherence 
length to the inter-carrier distance as a criterium of the 
BCS-Bose liquid crossover. The correct criterium is given 
by Eq.(5). 

II. SYMMETRY OF THE ORDER PARAMETER 



I S> £(0), so that the parameter- free expression for ep, 
Eq.(3), is perfectly applicable. 

A parameter-free estimate of the Fermi energy ob- 
tained using Eq.(3) is presented in the Table. The renor- 
malised Fermi energy of about 30 cuprates is less than 
100 meV, in particular if the degeneracy g > 2 is taken 
into account. That should be compared with the char- 
acteristic phonon frequency, which is estimated as the 
plasma frequency of oxygen ions, 

UirZ 2 e 2 N 

^ = y^4— (4) 

One obtains ujQ=SAmeV with Z = 2, N — 6/V ce u, 
M = 16 a.u. for Y Bci^CuzOq. Here V ce u is the volume 
of the (chemical) unit cell. The estimate agrees with the 
measured phonon spectra. As established experimentally 
in cuprates, the high-frequency phonons are strongly cou- 
pled with carriers. Therefore the low Fermi energy is a 
serious problem for the BCS (or Migdal-Eliashbcrg) ap- 
proach. Since the Fermi energy is small and phonon fre- 
quencies are high, the Coulomb pseudopotential u* is of 
the order of the bare Coulomb repulsion, fx* ~ fi c ~ 1 
because the Tolmachev logarithm is ineffective. Hence, 



The anomalous Bogoliubov average 

F. a /(r 1> r 2 ) = <(*.(r 1 )*y(r 2 ))) > 

is the superconducting order parameter both in the weak 
and strong-coupling regimes ( 4 , s (r) annihilates a carrier 
with spin s and coordinate r ). F ss '(ri,r2) depends on 
the relative coordinate p = r 1 — Y2 of two electrons of the 
pair, and on the center-of-mass coordinate R = (ri + 
V2j/2. Hence, its Fourier transform, /(k, K), depends 
on the relative momentum k and on the center-of-mass 
momentum K. In the BCS theory, where K =0 (in a 
homogeneous superconductor), the Fourier transform of 
the order parameter is proportional to the gap in the 
quasiparticle excitation spectrum, /(k, K) oc Ak. Hence 
the symmetry of the order parameter and the symmetry 
of the gap are the same in the weak-coupling regime. 
Under the rotation of the coordinate system, Ak changes 
its sign, if the Cooper pairing is d-wave. In this case the 
BCS quasiparticle spectrum is gapless. 

In the bipolaron theory the symmetry of the Bose- 
Einstein condensate is not necessarily the same as the 
"internal" symmetry of a pair |2^|. While the latter de- 
scribes transformation of /(k, K) with respect to the ro- 
tation of k, the former ("external") symmetry is related 
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symmetries, if 



n s [cos(7rm x ) ± cos^m^)] 



(8) 



Two order parameters, Eq.(8), are physically identi- 
cal because they are related by the translation transfor- 
mation, tpf (m x ,m y ) = ipj {m x ,m y + 1). Both have a 
d-wave symmetry changing sign, when the lattice is ro- 
tated by 7r/2, Fig.l. The d-wave symmetry is entirely due 
to the bipolaron energy dispersion with four minima at 
K^O. When the bipolaron spectrum is not degenerate 
and its minimum is located at T point of the Brillouin 
zone, the condensate wave function is s-wave with re- 
spect to the center-of-mass coordinate. The symmetry 
of the normal state (pseudo)gap has little to do with the 
symmetry of the order parameter in the strong-coupling 
regime. The one-particle pseudogap is half of the bipo- 
laron binding energy A/2, and does not depend on any 
momentum in zero order of the polaron bandwidth, i.e. 
it has an "s"-wave symmetry. In fact, due to a finite dis- 
persion of polaron and bipolaron bands, the pseudogap is 
an anysotropic s-wave. A multi-band electron structure 
can include bands weakly coupled with phonons which 
could overlap with the bipolaronic band. In this case 
CBG coexists with the Fermi gas, like A He bosons co- 
exist with 3 He fermions in the mixture of Helium-4 and 
Helium-3 01 • The normal state one-particle excitation 
spectrum of such mixtures is gapless. 



FIG. 1: d-wave condensate wave function, in the Wannier 
representation. The order parameter has different signs in 
the shaded cells and is zero in the blank cells. 



to the rotation of K. Therefore it depends on the bipo- 
laron band dispersion, but not on the symmetry of the 
bound state. 

As an example, let us consider a tight-binding bipo- 
laron spectrum comprising two bands on a square lattice 
with the period a = 1, 

= tcos(K x ) -t' cos(K y ), (6) 
E k = ~t'cos(K x ) +tcos(K y ). 

They transform into one another under 7r/2 rotation. 
If t, t' > 0, "x" bipolaron band has its minima at 
K=(±7r,0) and y-band at K = (0, ±7r). These four 
states are degenerate, so that the condensate wave func- 
tion ip s (m) in the site space, m = (m x , m y ), is given by 

^(m)=iV- 1 /2 J2 ke" Km . (7) 

K=(±7r,0),(0,±7r) 

where &k = are c-numbers at T = 0. The super- 

position, Eq.(7), respects the time-reversal and parity 



III. CHECKERBOARD MODULATIONS OF 
THE SUPERCONDUCTING GAP 

Independent observations of normal state pseudogaps 
in a number of magnetic and kinetic measurements, and 
unusual critical phenomena tell us that many cuprates 
may not be BCS superconductors. Indeed their super- 
conducting state is as anomalous as the normal one. In 
particular, there is strong evidence for a d-like order pa- 
rameter (changing sign when the CuO-i plane is rotated 
by 7r/2) in cuprates |2J|. A number of phase-sensitive 
experiments [25j provide unambiguous evidence in this 
direction; furthermore, the low temperature magnetic 
penetration depth [2(| |2tJ was found to be linear in a 
few cuprates as expected for a d-wave BCS supercon- 
ductor. However, different tunnelling spectroscopies, in 
particular c-axis Josefson tunnelli ng I28W , and some high- 
precision magnetic measurements 29] show a more usual 
s-like symmetry. Other studies even reveal an upturn 
in the temperature dependence of the penetration depth 
below some characteristic temperature (3^. Also both 
angle-resolved photoemission (ARPES) and tunnelling 
spectroscopies often show a very large energy gap with 
2A/T C ratio well above that expected in any-coupling 
BCS theory. 

Strong deviations from the Fermi/BCS-liquid be- 
haviour are suggestive of a new electronic state in 
cuprates, which could be a charged Bose liquid of bipo- 
larons. In the bipolaron theory the symmetry of the 
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Bose-Einstein condensate on a lattice should be dis- 
tinguished from the 'internal' symmetry of a single- 
bipolaron wave function, and from the symmetry of a 
single-particle excitation gap. As described above the 
Bose-Eistein condensate of bipolarons could be d-wave, 
if bipolaron bands have their minima at finite K in the 
center-of-mass Brillouin zone. At the same time the 
single-particle excitation spectrum is an anysotropic s- 
wave providing a possible explanation of conflicting ex- 
perimental observations. 

The two-dimensional patterns, Fig.l, are oriented 
along the diagonals, i.e. the d-wave bipolaron condensate 
is 'striped'. Hence there is a fundamental connection be- 
tween stripes detected by different techniques 0, Is2| in 
cuprates and the symmetry of the order parameter [2.1| . 
Originally antiferromagnetic interactions were thought to 
give rise to spin and charge segregation (stripes) [33| . 
However the role of long-range Coulomb and Frohlich in- 
teractions was not properly addressed. We showed that 
the Frohlich electron-phonon interaction combined with 
the direct Coulomb repulsion does not lead to charge seg- 
regation like strings or stripes in doped insulators, and 
the antiferromagnetic exchange interaction is not suffi- 
cient to produce long stripes either [34j . 

If (bi)polaronic carriers in many cuprates are in a liq- 
uid state, one can pose a key question of how one can 
see stripes at all. Actually, the bipolaron condensate is 
striped owing to the bipolaron energy band dispersion, 
as discussed above. In this scenario the hole density, 
which is about twice of the condensate density at low 
temperatures, is striped, with the characteristic period of 
stripes determined by inverse wave vectors correspond- 
ing to bipolaron band-minima. Such interpretation of 
stripes is consistent with the inelastic neutron scatter- 
ing in Y Ba2Cu^,Oj^s , where the incommensurate peaks 
were observed only in the superconducting state [35|. 
The vanishing at T c of the incommensurate peaks is in- 
consistent with any other stripe picture, where a char- 
acteristic distance needs to be observed in the normal 
state as well. On the contrary, with the d-wave striped 
Bose-Einstein condensate the incommensurate neutron 
peaks should disappear above T c , as observed. Also re- 
cent X-ray spectroscopy |3(| did not find any bulk charge 
segregation in the normal state of a high- Tc LaCuO^s 
suggesting an absence of in-plane carrier density modu- 



lations in this material above T c . 

The d-wave condensate reveals itself as a checkerboard 
real-space modulation of the hole density, Fig.l. Then 
the superconducting gap should be modulated as well. 
Indeed, as shown in Ref.j3g], the single-particle excita- 
tion spectrum of the bipolaronic superconductor can be 
obtained using the BCS excitation spectrum but with a 
negative chemical potential, 

e(k) - [(fc 2 /2m* + A p ) 2 + A C (T) 2 ] 1/2 , (9) 

where A C (T) 2 = constants \ 4> s \ 2 is a coherent component 
of the gap. This spectrum is quite different from the BCS 
quasiparticles because the chemical potential is negative 
with respect to the bottom of the single-particle band, 
(i = — A p . A single particle gap, A/ 2, defined as the 
minimum of e(k), is given by 

A/2 = [A 2 + A C (T) 2 ] 1/2 . (10) 
It varies with temperature from A(0)/2 = 

1 /2 

[A 2 + A c (0) 2 ] at zero temperature down to the 
temperature independent A p above T c . The spectrum, 
Eq.(10) was obtained for a homogeneous condensate. 
However, in the framework of a semiclassical approx- 
imation, we can apply Eq.(10) also with a modulated 
ip s (m), if the period of modulations is large enough 
|38| . As a result, the superconducting gap is modulated. 
Importantly, the checkerboard modulation of the gap 
have been observed in the tunnelling experiments with 
391 and without 

[HIHI applied magnetic field in Bi 

cuprates. 

In conclusion, I have shown that there is a real-space 
pairing in many cuprates due to the polaron collapse 
of the Fermi energy. Remarkably, the bipolaron theory 
provides an explanation of the d-wave order parameter, 
the anysotropic s-wave single-particle gap, and its stripe 
checkerboard modulations below T c within a single mi- 
croscopic approach. 
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